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Abstract 

We present a study of the parton distributions in transverse position or impact parameter space 
using a recent parametrization of u and d quark generalized parton distributions (GPDs) H{x, t) and 
E(x, t) at zero skewness. We make a comparative study between different parametrizations and discuss 
the region of vahdity of the positivity condition. 
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Introduction 

Generalized parton distributions (GPDs) contain a wealth of information about the nucleon 
structure (see [jj for example). Unlike the ordinary parton distributions (pdfs) which at a 
given scale depend only on the longitudinal momentum fraction x of the parton, GPDs are 
functions of three variables, x, ( and t where the so-called skewness ( gives the longitudinal 
momentum transfer and —t is the square of the momentum transfer in the process. These 
are called off-forward parton distributions. The GPDs give interesting information about the 
spin and orbital angular momentum of the quarks and gluons in the nucleon as well as their 
spatial distribution. They are experimentally accessed through the overlap of deeply virtual 
Compton scattering (DVCS) and Bethe-Heitler (BH) process as well as exclusive vector meson 
production [11]. Data has been obtained at HERA collider, by the HI 0, 3 and ZEUS Q, 3 
collaborations and HERMES Q] fixed target experiment. DVCS experiments are also being 



done at JLAB 
muon beams 



Hall A and B 



7|. COMPASS at CERN has programs to access GPDs through 
. Experimental observables, however involve a convolution of GPDs, and so 
modeling GPDs is interesting. GPDs reduce to ordinary forward parton distributions (pdfs) in 
the forward limit. Moments over x give nucleon form factors. These act as useful constraints 
to model the GPDs. A large number of models or parametrizations have been proposed for 
GPDs. Here we do not plan to review all of them but mention only those that are relevant for 
us. A detailed list of the main lines of approach and their present status with respect to the 
data can be found in 3. Moments of GPDs have been calculated on lattice as well. In 10], 
GPDs at zero skewness ( have been parametrized in a Regge type model at small x and using 
a modified Regge parametrization at large x and t. It was found to describe the basic features 
of proton and neutron electromagnetic form factors. Similar small x Regge type behaviour was 
used in the modeling of pion GPDs in Another parametrization of GPDs at zero skewness 
is given in [l^. Here a Regge motivated x dependence at small x was interpolated to large x 
region. The GPDs at the input scale were fitted to the experimental data on Dirac and Pauli 
form factors. An exponential t dependence was used. The form of the t dependence was found 
to be unchanged by the scale evolution. In [13] a recent parametrization was proposed for 



zero skewness and an extension to non-zero ( was done in 



ij]. At the input scale the GPDs 



were parametrized by a spectator model term multiplied by a Regge motivated term. The 
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15, 



lej. A recent fit to the DVCS data 



parameters were obtained by fitting the forward pdfs and form factors. 

For non-zero (, the GPDs have to satisfy an additional constraint, namely polynomiality. 
In certain models, for example using the overlap of light-front wave functions (LFWFs), it is 
very difficult to obtain a suitable parametrization of the higher Fock components of the wave 
function in order to get the polynomiality of GPDs. Polynomiality is satisfied by construction 
only if one considers the LFWFs of simple spin 1/2 objects like a dressed quark or a dressed 
electron in perturbation theory instead of the proton 
at small Bjorken x from HI and ZEUS was done in 
representation of the DVCS amplitude. However, to get the GPDs one has to do an inverse 
Mellin transform, and a knowledge of all moments are required for that. 

At zero skewness (, if one performs a Fourier transform (FT) of the GPDs with respect to 
the momentum transfer in the transverse direction A^, one gets the so called impact parameter 
dependent parton distributions (ipdpdfs), which gives how the partons of a given longitudinal 
momentum are distributed in transverse position (or impact parameter b±) space. These obey 
certain positivity constraints and unlike the GPDs themselves, have probabilistic in terp retation 



17|, using the conformal Mellin-Barnes 



19|. Due to 



[I8l |. These give an interesting interpretation of Ji's angular momentum sum rule 
rotational invariance, the same relation should hold for all components of the angular momen- 
tum Jq. In the impact parameter space, the relation for has a simple partonic interpretation 



for transversely polarized state 



201 1 : the term containing E{x,0,0) arises due to a transverse 



deformation of the GPDs in the center of momentum frame. The term containing H{x, 0, 0) 
is an overall transverse shift when going from the transversely polarized state in instant form 
(rest frame) to the front form (infinite momentum frame). On the other hand, in 2l|], real 
and imaginary parts of the DVCS amplitudes are expressed in longitudinal position space by 
introducing a longitudinal impact parameter a conjugate to the skewness (, and it was shown 
that the DVCS amplitude show certain diffraction pattern in the longitudinal position space. 
Since Lorentz boosts are kinematical in the front form, the correlation determined in the three- 
dimensional b±,a space is frame-independent. As GPDs depend on a sharp x, the Heisenberg 
uncertainty relation restricts the longitudinal position space interpretation of GPDs themselves. 



It has, however, been shown in 



22l | that one can define a quantum mechanical Wigner distri- 



bution for the relativistic quarks and gluons inside the proton. Integrating over k and k- 
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one obtains a four dimensional quantum distribution which is a function of r and where r 
is the quark phase space position defined in the rest frame of the proton. These distributions 
are related to the FT of GPDs in the same frame. This gives a 3D position space picture of 
the GPDs and of the proton. 

In this series of work, we plan to investigate the GPDs for the proton in transverse and 
longitudinal position space. In this first paper, we study the recently parametrized form in [l^ 
in impact parameter space and make a comparative study with other models. 

Parametrization of the GPDs 

We consider the parametrization in [l^ for the GPDs : 
Set I 

H\x,t) = Gj;;.(a;,t)x-"'-^^(i-^)'''* (1) 
E\x,t) = KG^,(a;,t)a;-"'-^2(i-")'^* (2) 



Set II 



All parameters except for pi and p2 are flavor dependent. The function G has the same form 
for both parametrizations, I and II: 



where 



D{x,\iA_) = k'^-m^, (6) 



e = xM' - -^Ml - ^ (7) 
1 — X 1 — X 

k'^ = xM^ Ml - ^ (8) 

l-x l-x ^ ' 
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and 

^(k\X) = ^^,. (9) 

Here the skewness variable is taken to be zero, in other words, momentum transfer is in 
the transverse direction, t is the invariant momentum transfer squared, t = — A^, and x is the 
fraction of the hght cone momentum carried by the active quark, k being its momentum. The 
mass parameters are m, the struck quark mass, and M, the proton mass. The normahzation 
factor includes the nucleon-quark-diquark coupling, and it is set to A/" = 1 GeV^. 

The u and d quark contributions to the anomalous magnetic moments are: 

, k"^ = -2.03, hi q = d 
k = k'^ = { . (10) 

k"/2 = 1.67/2, for g = M 



The parameters are listed in [13] for both the sets. The parameters M^, A'' and a'', q = u,d, 
obtained at an initial scale Ql {Q^ = 0.094 GeV^), and they are the same for both Sets I and 
II, in Set I they are by definition the same for the functions H and E (see Eqs. (|T|2l) ). The 
parameters in this model are obtained by fitting the experimental data on the nucleon electric 



and magnetic form factors (see 13| for reference). For the forward limit, Alekhin 23| leading 
order (LO) pdf sets were fitted within the range 10-^ < X < 0.8 and 4 < < 240 GeV^ by 
valence distribution and using the baryon number and momentum sum rules. The input scale 
Ql = 0.094 GeV^ is obtained parameter in this model. The low value of Qq results from 
the requirement that only valence quarks contribute in the momentum sum rule. 

The above phenomenologically motivated parametrization of the GPDs H{x,t) and E{x,t) 
at zero skewness ( was done using a spectator model calculation at the low input scale. The 
spectator model has been used for its simplicity and for the fact that it is flexible enough 
to predict the main features of a number of distribution and fragmentation functions in the 
intermediate and large x region. The spectator mass is chosen to be different for different quark 
flavor GPDs. However, similar to the case of pdfs, the spectator model is not able to reproduce 
quantitatively the small x behaviour of the GPDs. So a 'Regge-type' term has been considered 
multiplying the spectator model function G^^^. The parameters were obtained by fitting the 
form factors and forward pdfs. Two versions of the parametrizations were used and are given 
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by set I and set II. The GPD E is unconstrained by the data on forward pdfs, so in set II an 
additional normalization condition has been imposed 

"1 

dxEq{x,t = 0) = k'^ (11) 



with the experimental values of and k,'^. 

Although if" and H'^ are similar in behaviour in both sets of parametrization, the main 

I > 

difference is in the behaviour of E'^ and E'^ (see Fig. 9 of [13]) at the input scale. E^ is a 
slowly increasing function of x for set I, and for set II, it increases rapidly as x approaches 
zero. Ed has a peak at larger value of x for set II. It was found that the difference between the 
sets decreases if one evolves the GPDs to higher values of (scale). It is to be noted that 



the parametrizations for [idl and [12 1 are at different input scales compared to 1131. However 
the GPDs are evolved to the respective scales of pi| and and a comparison is provided in 
Figs. (13-16) of [l^. For and H'^ although there is agreement for low — t, for higher —t the 
results differ qualitatively and quantitatively. For i?" and E'^ the disagreement is also at lower 
-t. 

Again, even 

at = 4 GeV^ there is significant difference between set I and set II. H'^ and 
H'^ agrees with lattice calculations of [24] at t = —0.3 GeV^, but the qualitative behaviour is 
different and even is outside the error band of lattice calculations for higher values of | t | . 

Parton distributions in impact parameter space 



Parton distribution in impact parameter space q{x,b) is defined as 181 ]: 



q{x,b) = ^ I d'Ae''^^-'^H{x,t) 

e(x,b) = At /" d^Ae-'^^-'^Eix^t). (12) 
An^ J 

These functions have the physical interpretation of measuring the probability to find a quark 
of longitudinal momentum fraction x at a transverse position b± in the nucleon. Here b = \ b± \ 
is the impact parameter which is the transverse distance between the struck parton and the 
center of momentum of the hadron. b is defined such a way that ^2- Xib-i = where the sum 
is over the number of partons. The relative distance between the struck parton and the 
spectator system provides an estimate of the size of the system as a whole. 
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In Fig. 1-2 we have plotted q{x,b) and e{x,b), both for u and d quarks and for set I. The 
values of the parameters used are from [13'] at the input scale. For small and medium x, e'^{x, b) 
is larger in magnitude than e"(x, b). The peak shifts to higher x as 6 decreases. This means that 
the d quark dominates in the proton helicity flip distribution. However, u quark contribution 
dominates in the helicity non-flip q{x,b). e'^ is negative whereas is positive, similar to the 
model in [l^. However, in the model we study, e'^{x,b) is comparable with or even larger in 
magnitude than q'^{x,b), unlike in 10|, where it is much smaller at the input scale. In the 



parametrization of [10], the t dependence is only in the argument of the exponential and the 
Fourier transform is simpler and can be obtained analyically. The resulting parton distributions 
in the impact parameter space are larger in magnitude than what we get in this work using 
the parametrization of [13] which may be due to the different scale Qq in the plots. In 121 ]. 
2 — D distributions in the b^ — by plane are plotted. The smearing of the quark distributions in 
the transverse impact parameter plane decreases as x increases, which means that the parton 
distributions are more localized for higher values of x. Similar behaviour is observed in the 
model of [l^. As x approaches 1, the transverse width of q{x,b) should vanish 18]. In this 
limit q{x, b) should have a very peaked transverse profile, as H{x^ t) is independent of t when 
X — > 1 as the active quark carries all the proton momentum no matter what t is. However, 
from Fig. 1 we see that as x ^ 1, the peak of the distribution decreases. 

g(x, b) both for u and d quarks in set II are the same as in set I. In Figs. 3, we have plotted 
e(x, b) for u and d quarks where the parameters are as in set II. e(x, b) has a different behaviour 
compared to set I. The peak of e" is shifted to very small value of x and e"(x, 6) decreases 
sharply as x increases. That means at larger x, d quark dominates in e(x, 6). 

The Fourier transform of the GPD E{x,t) plays an important role when instead of an 
unpolarized target we have a transversely polarized target. In other words, it has a probability 
interpretation in the transversity basis rather than the helicity basis. For a state polarized in 
the X direction, parton distribution in the impact parameter space becomes [isl 

g (x,&) = ^(x,fe)-^^^ (13) 

This means that the GPD E{x,t) causes a transverse shift of the quark distribution in a 
transversely polarized target. For a state polarized in x direction the shift is in the y direction 
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and so on. The magnitude of the shift is given by \ ^g^'^-* |. The average displacement of 
the shift is given by 

_ jd%hyq^{x,h) _ 1 E^(x,0) 
^ ' ^ jd%q^{x,b) 2MHi{x,Qy ^ ' 

The distance between the struck quark and the spectator system is given by 

For d quarks, s'^(x) is larger in magnitude than u jl2i|. The transverse shift depends on the set 
of parameters used in the model considered here [l3|. 

The parton distributions in impact parameter space should obey the positivity condition 



25| 



This follows from the fact that q^{x, b) which is the unpolarized parton distribution for trans- 
versely polarized proton should have probability interpretation. 

The positivity relation in effect puts an upper bound on b. In the model of lOj, it was found 
that for u quarks, the positivity bound was satisfied over most of the x region, considering 
that the GPDs are vanishingly small for values of b larger than the nucleon size. For d quarks, 
a violation of positivity was observed which is more pronounced at larger values of x and b. 



However the violation is small. In the model of 



12], the positivity condition was used to 



constrain the behaviour of E{x,t = 0) at large values of x. In Fig. 4 we have plotted q{x,b) 



as well as \ ^^^^ | for u and d quarks for both the parametrizations in [ISl], as functions of 
b. We see that in set I for u quarks, positivity is mildly violated for very large values of x and 
for b > 0.5 GeV^^. For u quarks in set II there is no violation. For d quarks in set I, positivity 
is not violated for any b values for a; < 0.1. It is violated at larger b values or comparatively 
larger x. Violation is more as x increases. For d quarks in set II, positivity is largely violated 
for large x, b values. To see the x dependence of the violation of positivity, in Fig. 5 we plot 
the same quantities as in Fig. 4, as functions of x, for fixed values of b. For u quarks in set 
I positivity is violated for b ~ 0.4 GeV~^ for large values of x. For higher b values, violation 
is even at smaller values of x. Violation of positivity is not seen in u quarks in set II. For d 



quraks in set I, there is larger violation of positivity and it starts already at x = 0.2 for large 
b. For d quarks in set II there is again a large violation of positivity that starts at a; ~ 0.3 for 

b = 5 GeV-\ 



Conclusion 

In this paper, we have studied the GPDs H{x,t) and E{x,t) for zero skewness in a recently 
parametrized form in transverse position or impact parameter space. We present a com- 
parative study between several models. A violation of the positivity condition was observed 
in certain range of x, b. This puts additional constraint on the kinematical region where this 
parametrization is to be used. This depends on the set of parameters used in the model. A 
new fit of the parameters with this additional constraint may improve the model. As extension 
to nonzero ( was proposed in [ij]. In a future work, we plan to investigate the GPDs with 
nonzero skewness both in transverse and longitudinal position spaces. In particular, a study 
in longitudinal position space is interesting to understand the origin of the observed diffraction 
pattern in the DVCS amplitude in a simple QED model j2l| . 
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FIG. 1: (Color online) Plots of (a) q^{x, b) vs b =\ b± \ for fixed values of x, (b) g"(x, b) vs x for fixed 
values of b, (c) same as in (a) but for q"^, (d) same as in (b) but for q'^. Parameters are as in set I. b 
is in GeV-^ 
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(c) b (d) X 



FIG. 2: (Color online) Plots of (a) e'^{x,b) vs b = \ b± \ for fixed values of x, (b) e"(x, b) vs x for fixed 
values of b, (c) same as in (a) but for -e'^, (d) same as in (b) but for —e'^. Parameters are as in set I. 
b is in GeV^^ 
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(c) b (d) b 

FIG. 4: (Color online) Plots of (a) q^{x,b) and y^{x,b) = \ q^'^^ \ vs. b for set I, (b) the same 
as in (a) but for set II, (c) q'^{x, b) and y^{x, b) = \ ^^'^^ \ vs. b for set I, (d) the same as in (c) 
but for set II. b is in GeV^^. 
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FIG. 5: (Color online) Plots of (a) q^{x,h) and y^{x,h) = \ ^ q^' ' \ vs. x for set I, (b) the same 
as in (a) but for set II, (c) q'^{x, b) and y'^{x, b) = \ g^'^-* | vs. x for set I, (d) the same as in (c) 
but for set II. b is in GeV^^. 
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